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ABSTRACT 

Sunyaev-Zeldovich (SZ) surveys find massive clusters of galaxies by measuring the inverse Compton 
scattering of cosmic microwave background off of intra-cluster gas. The cluster selection function from 
such surveys is expected to be nearly independent of redshift and cluster astrophysics. In this work, 
we estimate the effect on the observed SZ signal of centrally-peaked gas density profiles (cool cores) 
and radio emission from the brightest cluster galaxy (BCG) by creating mock observations of a sample 
of clusters that span the observed range of classical cooling rates and radio luminosities. For each 
cluster, we make simulated SZ observations by the South Pole Telescope and characterize the cluster 
selection function, but note that our results are broadly applicable to other SZ surveys. We find that 
the inclusion of a cool core can cause a change in the measured SPT significance of a cluster between 
0.01%-10% at z > 0.3, increasing with cuspiness of the cool core and angular size on the sky of the 
cluster (i.e., decreasing redshift, increasing mass). We provide quantitative estimates of the bias in 
the SZ signal as a function of a gas density cuspiness parameter, redshift, mass, and the 1.4 GHz radio 
luminosity of the central AGN. Based on this work, we estimate that, for the Phoenix cluster (one 
of the strongest cool cores known), the presence of a cool core is biasing the SZ significance high by 
^6%. The ubiquity of radio galaxies at the centers of cool core clusters will offset the cool core bias 
to varying degrees. 

Subject headings: cosmology: observations - galaxies: cooling flows - galaxies: clusters 


1. INTRODUCTION 

Galaxy clusters are potentially powerful t ools to study 
dark energy and cosmology (see review by Allen et al. 
20111; whereas standard candles test cosmology on ho- 


mogeneous scales, clusters trace the growth of inhomo¬ 
geneity via N{M,z), the cluster number density at a 
given mass and redshift. Accurate cluster cosmology re¬ 
quires not only reliable estimators of clusters properties 
such as total mass (i.e., M^qq), but also an accurate es¬ 
timation of the survey selection function. 

A new technique of assembling a nearly unbiased sam¬ 
ple of galaxy clusters uses the Su nyaev-Zerdovich (SZ) 
effect ( Sunyaev fc Zeldovich|1972' ): the inverse Compton 
scattering of photons as they pass through the hot intr¬ 
acluster medium, leading to distorted-spectrum patches 
in the cosmic microwave background (CMB). One attrac¬ 
tive characteristic of SZ surveys is that the SZ signal is 
virtually redshift independent, as opposed to the ^ l/d| 
flux dimming in infrared, optical, and X-ray surveys. Re¬ 
cently, large samples (>100) of galaxy clusters have been 
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Reichard 
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many of the most massive known galaxy clusters (e.g 
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Many potential systematics have been considered in 
SZ-selected cluster surveys. Simulations have shown that 
SZ surveys are expected to be relati vely insensiti ve to the 
effects of no n-gravitational p hysics (|Nagai|2006|, projec¬ 
tion effects ( Shaw et al.|2008|), a nd the dynamical state of 
the cluster ( Krause et al. [21)12 1. The presence of a “cool 


core” — a central region of over-dense gas accompanied 
by a drop in temperature — could also conceivably bias 
mass estimates. Finally, contamination from radio-loud 
active galactic nuclei (AGN), generally located in the 
brightest cluster galax y, could introduce bias by “filling 
in” the SZ decrement (Sayers et al.||2012|). 

Gool core bias is also of particular astrophysical inter¬ 
est, as estimates of the cool core fraction as a function of 
redshift might lead to valu able insights into the cooling 
flow problem (see review by Fabian 1994). Recent studies 
have found evidence for a lack ot dense, cool c ores in the 
centers of high-redshift galaxy clusters (e.g., |Vikhlinin 


et al. 1120071 ISantos et al.||?010| |McDonaId||2Ull[ [McDon¬ 


ald et al. 2013|). These results may be evidence for an 
evolution in the cooling/feedback balance, or may be a 
result of biases in X-ray and optical surveys (e.g., Sem- 
ler et al. 2012). Thus, in order to fully understand the 
evolution ot cooling and heating processes in the cores 
of galaxy clusters, we must understand how our sample 
selection affects the observed cool core fraction. 

Throughout this paper, we take the fractional mat¬ 
ter density Dm = 0.27 and the fractional vacuum den¬ 
sity Da = 0.73, and neglect radiation D^ and curvature 
Df.. We introduce two dimensionless forms of the Hub¬ 
ble parameter, h = iJo/(100kms“^Mpc“^) = 0.7 and 
E{z) = H{z)/Ho. We define the cluster radius rsoo such 
that the average matter density interior to r^oo is 500 
times the critical density pcr of the universe as well as 
the cluster mass M^qo = |’’■7■5oo ^ SOOpcr- 












































2 


2. METHODS 

In order to measure the influence of cool cores and 
radio-loud AGN on the SZ signal, we generate one¬ 
dimensional (spherically-symmetric) mock galaxy clus¬ 
ters with properties drawn from existing samples of high- 
mass, relaxed clusters, and add a perturbation either to 
the gas density (cool core) or mm flux (radio source). 
Simulated observations of these mock clusters with the 
South Pole Telescope are generated, including realistic 
noise and background, yielding a realistic SZ signal-to- 
noise measurement. 


2.1. Constructing Mock Galaxy Clusters 

Since we are mostly interested in the SZ signal of our 
mock clusters, the key ingredient in our cluster models is 
the pressure profile P{r). However, we must also model 
the cluster density profiles p{r) since we wish to calculate 
the change in the pressure profile 5P due to a density per¬ 
turbation 5/0, which in general is not 5P oc 5p as the tem¬ 
perature profile also shifts T ^ T — ST. Consequently, 
we need an additional constraint on our clusters to cal¬ 
culate SP, which we take to be hydrostatic equilibrium: 
^ = —p^, where $ is gravitational potential. Since 
there is no evidence that the underlying dark matter dis¬ 
tribution i s affected by the bary onic processes driving 
cool cores ( {Blanchard et al.||2013 ), 5<i) is straightforward 
to calculate, and using P{r 5 oo) as a boundary value, it is 
possible to calculate P + SP. 

The details of our density and pressure parameteriza- 
tions are now presented. For each non-cool core cluster, 
we draw a cluster mass M^qq and a redshlft from uniform 
distributions (see Table [^, with ranges motivated by the 
observed r anges in the SPT 2500 deg^ survey (Bleem 
et al.|[MT4 ). The gas density profiles were then modeled 
by a slightly modified version of the functional form of 
Vikhlinln et al. (2006): 


TlpTlQ — 


nl{r/rc + S)~ 


[1 + (r/r,)2]3/5-«/2 


( 1 ) 


where the various length scales and slopes were sampled 
from realistic ranges for massive galaxy clusters by a 
Monte Carlo process, which we detail in Table For 
this work, we have “r egulated” the parameterization of 
Vikhlinin et al. (2006) by inserting a small factor S = 0.1 
to obtain a finite density at r = 0. This is necessary 
as we will eventually extract an SZ signal from the en¬ 
tire cluster, whereas the behavior of the profile very near 
r = 0 is not of interest to observational works due to the 
limited spatial resolution of any survey. 

For a pure non-cool core, a = 0 and equation ([^ re¬ 
duces to the simpler form 


rin 


[1 -b {r/rc)^f^ [1 -b {r/rs)^Y^^ 


( 2 ) 


where a, (3, and e are slope parameters for the r r^, 
'^c ^ f fs and r ^ rs radial regimes. We convert UpUe 
to the gas mass density via the relation p = mpUeAlZ 
where Z = 1.199 is the average nuclear charge, which is 
greater than unity because of ionized elements heavier 
than hydrogen, A = 1.397 is the average nuclear mass, 
which is greater than Z due to the presence of neutrons, 
and Uf, = Zup is the electron density in terms of the 


positive nuclei density Up. The density normalization 
no was chosen to enforce an average gas fraction /gas = 
0.125 within rsoo. 

For the unperturbed, non-cool core clusters, we adopt 
the simple temperature profile 


T{r) 

Tn 


= 1.35 


1-b 


O.brsoo 


-0.45 


(3) 


which is the “universal” temperature profile of |Vikhlinin| 
et al. (2006) with the term accounting for the central 
temperature drop suppressed. 

For each non-cool core cluster, we generate 14 progres¬ 
sively cuspier cool-core clusters. We start by duplicat¬ 
ing each non-cool core cluster and shifting p —>■ p -b Sp. 
To avoid adding more free parameters (a cooling radius, 
etc.) which may be correlated with the other model pa¬ 
rameters, we simply increment a to steepen the density 
profile towards the center of the cluster. Each successive 
cluster has an = ctn-i -bO.25. The density normalization 
no was also reduced by an appropriate factor to keep the 
gas mass within r 2 oo constant. 

Finally, the assumption of hydrostatic equilibrium was 
used to recalculate the pressure profile of the cool core 
cluster P —)■ P -b SP: 


p + SpAP {p + Sp)GSM 

p df p2 ^ 

We vary P(r 2 oo) until the “internal energy” oc / PdV of 
the cool core cluster out to r 2 oo agrees with its non-cool 
core counterpart. Note that in general, this means that 
SP{r 2 oo) ^ 0, though at r 2 oo, SP < P. 

The above methods were used to generate 200 x 15 
mock clusters which ranged from pure non-cool core 
a = 0 to strong cool core a = 3.5. We chose the maxi¬ 
mum value a = 3.5 as even strong cool core clusters like 
Phoenix have a < 3.5. 


P + SP = / df 


2.2. Model Validation 

Additional quantities were derived to demonstrate how 
the wide range of observed cluster properties are cap¬ 
tured by our methodology. 


1. Temperature. Our derive d temperature profiles are 
displayed in Figure [I(b)| The cool core profiles ex¬ 
hibit a central temperature drop similar to the av¬ 
erage co ol core profiles observ ed in X-ray selected 
samples (Vikhlinin et al.||2006). The characteristic 
temperatures of our clusters range from ~ 1 keV 
to ~ 30 keV, which is in agreement with observa¬ 
tions of massive clusters. Furthermore, Tmin/To, 
which corresponds to the ratio of the central cool 
core core temperature to the a = 0 ce ntral tem¬ 


perature, ra nge from 0.1 to 0.85 in the (Vikhlinin 


et al. 2006) clusters, which is completely covered 


by our simulated clusters. Similarly, our simula¬ 
tions cove r the range of temper ature profiles ob¬ 
served in (McDonald et al. 2014), which find that 
the fractional temperature in the inner O.OdPsoo is 
Tcore/Tsoo = 0.74^09o4 for high redshift, cool core 
clusters and higher for low redshift cool core clus¬ 
ters. 
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TABLE 1 

Parameters. All radii are in units of kpc. We report la uncertainties. 


Parameter Range 


Distribution Notes 


Afsoo/Al© [1 X 10^^, 2.1 X 10^®] Log-Uniform 

z [0,2] Uniform 

SPT 2500 deg2; 

31eem et al. 

(2014 

) 

) 

SPT 2500 deg2; 

Bieem et al. 

(2014 

/3 (0.28,0.88,0.71) Triangle^ 

e 7.1-6/3'’ — 

^c/^’soo (0.067,0.26,0.47) Triangle® 

logr^/rsoo 1.21 — 

McDonald et al. 

(2013 

Vikhlinin et al. (2006) 

Andersson et ai. 

(2011 


iVlcUonald et al. 

(2013 

McDonald et al. 

(2013 


long Gaussian tails that lead to unphysical clusters. We report the (min, max, mode) ot the distribution in the 
range column. 

* The mock SZ signal is virtually e-independe nt. We choose e to let dlogP/dlogr = —0.90 at large radii, as 
given by the universal temperature profile of|Arnaud et al.1(|2010[. 






Fig. 1. — Clockwise from the upper left: density, temperature {Tx = f n^T AV j f n^dV), cooling time, and entropy as a function of 
cluster radius for a typical mock cluster. All of the displayed profiles were obtained by cloning a single non-cool core core cluster and 
adding density perturbations. Thes e plots represent 1/200 th of the total dataset. In the upper right plot, the thick dashed line represent 
the universal temperature profile ofIVikhlinin et al.N2006ll. In the lower left plot, the thick dashed line represents the predicted entropy 
profile slope from non-radiative simulations otfvoit et ah] ([2005^. 


2. Entropy Parameter. The entropy profiles of our 
mock clusters are derived in an effort to show their 
resemblance to observed clusters. Following the X- 
ray survey conventions, an entropy parameter K = 

-2/3 


ksTue 


is introduced. 


We pl ot thi s parameter 
Consistent 




with 

ICavagnolo et al. 

(2009) 

Hudson et al. 

(20101 3 

and 

McDonald et al. 

2013), the cool core clusters 


have a central entropy Kq < 100 keV cm^, whereas 
for non-cool core clusters Kq > lOOkeVcm^. The 


slope of our entropy profiles outside of the core 
is also in reaso nable agreement with the accretion 
shock model of Voit et al. (2005), which represents 
the “baseline” entropy profile ot clusters if radiative 
and other non-gravitational processes are ignored. 

The coo l ing ti me tcooi has been 


shown by Hudson et al. (2010) to be a quantity 
that can be used to se grega te cool cores from non¬ 


cool cores. In Figure 1(d) tcooi = 3fcT/ (ngA) is 
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plotted as a function of rad ius. Here, A(T, Z) is 
the cg oling function given by [Sutherland fc Dopita 
(19931. Most non-cool cores have central cooling 
times tcooi.o ^ ^/Hq whereas the cool core clusters 
have cooling t imes frnni n ~ I Gyr, consistent with 


the findings of |Hudson et al. (2010). 


2.3. Sunyaev Zel’dovich Maps 

For each mock cluster, we constructed an SZ map with 
an angular resolution of AO = 0.125 arcmin, which was 
converted to spatial resolution Ar 2 D = dA{z)A9 where 
dA^z) is the angular diameter distance to the cluster. 
The SZ line of sight integral was then calculated at each 
pixel: 


AT/Tcmb = ^^dl IT) 


atneksT I 


( 5 ) 


m-eC^ j 

/sz(i^,T) = [l + 6,,iX,T)] (6) 

where X = hv/ksTcMB, is the Thomson cross section, 
TOe is the mass of the electron, c is the speed o f light, and 
is the relativistic correction as given by jltoh et ah] 
(|1998|). For our purposes, it sufficed to evaluate dsz{T) at 
We made simulated SZ maps at 97.6 and 152.9 
GHz, corresponding to the effective freq uency of the SPT 


obser ving bands used for cluster-finding (Carlstrom et al. 
2011 1 . Our mock SZ maps were generated at a resolution 
<C 1' and then later convolved with the SPT beam which 
has an effective FWHM of ^ 1.2'. 

The mock SZ maps were then spatially hltered in a way 
identical to the SPT cluster-finding algorithm used in 


Reichardt et al. J 2013). A multi-frequency matched spa¬ 
tial hlter ( [Haehnelt I'egmark 1996 Melin et al.||2006 | 
was applied to the SZ maps in F'ourier space, which ac- 
counted for the cluster gas profile, and the other sources 
of astrophysical signal and instrumental noise expected 
in the SPT maps. The spatially hltered signal at the 
true cluster position was then re-normalized so that it 
w as equivalent to the unbia sed SPT signihcance, used 


Vanderlinde et al. (2010), which would effectively cor- 
respond to the signal-to-noise of cluster in a SPT map. 

These SZ maps were run through the SPT cluster- 
hnding pipeline, which adds noise in the form of point 
sources and CMB anisotropy maps, before attempting to 
detect clusters. For each cluster, the SPT pipeline then 
reported C , the “unbiased sign ihcance” (see Vanderlinde 


et al.|2010 Benson et al.|2011 ) that roughly corresponds 
to the cluster signal-to-noise ratio. 

2.4. Active Galactic Nuclei 

In order to test the effects of including radio-loud AGN, 
we created mock clusters with the following grid of pa¬ 
rameters: 

• M 500/2 X 10^‘^Mq = 10°, 10^/2, IQi 

• z = 0.3,0.4,0.5,0.7,1.1,1.4,1.7 

and we include a radio-loud AGN at the cluster center, 
with a range of luminosities: 

• logLi .4 = 23 - 29, AlogLi ,4 = 1.0 


where Li ,4 is the 1.4 GHz luminosity in WHz“^. All 
other parameters are set to their median values. The 
range in log Li ,4 covers the most luminous brightest clus¬ 
ter ga laxy (BGG) in the 152 X-ray BCG sample of Sun 
(20091, allowing for unusually shallow spectrum systems 
that would generate a large bias. 

We convert these luminosities to the SPT observing 
frequencies by assuming a spectral slope of as = 0.89, 
where hux scales with . This slope was the median 
spectrum between 1.4 and 30 GHz for radi o galaxies in 
a samp le of 45 massive clusters observed by [Sayers et al.j 
( 2012 ), however it is always possible to reinterpret the 
results for different choices slopes. We discuss deviations 
from this assumption in Section 3.1. As a control, a point 
source is not added to one cool-core/non-cool core cluster 
at each mass and redshift. 

3. RESULTS 


In Figure we show the 
SPT detection significance 
core cuspiness (Vikhlinin et 


(0500 = r 5 oo/dA[z}), and 


fractional change in the 
as a function of the cool 
angular size 
< 3', the 


al. 
redshift.' 


2007[), 

At 0500 


SPT detection significance changes by AC/C = (Ccc ~ 
Cncc)/Cncc ^ 0-1 for any type of cool core at any red- 
shift. The cool core bias is most pronounced for low- 
redshift, high mass systems, where the angular size of 
the SPT beam is much smaller than the angular size of 
the cluster so that the inner pressure profile becomes im¬ 
portant. In principle, the resulting redshift-dependence 
in the bias might lead to a ~ 2% bias towards cool core 
clusters from z ^ 1.5 to z ~ 0.3 in the observed cool 
core fraction CCF = Ncooi core/ACiusters evolution, but 
these effects should be negligible when comp ared to the 
^ 30% evol ution in CCF reported in e.g. 

eraIl([M3). 


McDonald 


In real observational studies it may be difhcult to de¬ 
termine the density shape parameter a, as this requires a 
multi-parameter fit to the UpUe prohle and assumptions 
on the shape of the cool core/non-cool core profiles. To 
ameliorate the situation, we calculate the “cuspiness” pa¬ 


rame ter 7 = —dlogp/dlogr|r'=o.o 4 r 5 nn (Vikhlinin et al. 
2007) of all simulated clusters. Since this quantity 
does not refer to the details of our parameterizations, 
it serves as a model-independent indicator of cool core 
strength; correlations between these quantities and A(/( 
are shown in Figure The bias is well-captured by the 
following relation: 


(<5C/C) = 10 


-2.08±0.01 


(Ay) 


0.85±0.01 


0500 (^) 


arcmin 


1.52±0.02 


( 7 ) 

where Ay = y — y|a=o « y — 0.1. This expression gives 
a practical method for estimating the SZ bias of a given 
cool core cluster, assuming that z > 0 . 1 . 

3.1. Radio bias 


Assuming a typical spectral slope of as = 0.89 (Say- 
ers et al.[[ 2012 [), only sources with logLi ,4 > 26 result in 
[AC/CI > 10% for z > 0.3 and M > 3 X lO^^. Conse¬ 
quently, none of the 152 radio-loud AGN listed in [Sun| 
(2009) would produce a bias in excess of 10% if located 


in massive clusters. More generally, the radio-loud point 
source bias can be readily estimated by the following fit 
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log 1O-2-O8 yO.8501.52 




Fig. 2.— Dependence of the SZ bias on the cuspiness (7 = ~ dlogr redshift (z), and angular size (©500 = '^ 50 o/^a{^)) for 3000 mock 
clusters. In the left panel, the point color corresponds to redshift of the mock cluster. In the central panel, we color by the angular size of 
the cluster. The clear separation into colored bands in this plot suggests that the bias is more fundamentally tied to the angular size than 
it is to the redshift. In the right-most panel, we show the edge-on projection of the best-fit three-dimensional plane (Equation 7). The 
~0.5 dex scatter is primarily due to variations in the cluster gas density profile, and the imperfect classification of cool core strength based 
on the 7 parameter. 








Fig. 3.— Dependency of A^/<^ on M^qq, redshift, radio luminosity, and spectral slope. The contours represent lines of constant log(A(^/(^). 
If an AGN exists in the shaded region, it will overpower the SZ signal. The large black dot represents a Phoenix-like cluster at multiple 
redshifts. Even at 2: = 0.3, a Phoenix-like cluster will exhibit only a 1% bias. The most extreme s ystems exhibit a ~ 10% bias at 2 ~ 0.3. 
Red dots represent a sam ple of X-ray selected massive ell iptical galaxies from |Dunn et al.| l|2010|l and blue dots represent extreme AGN 
hosted in massive clusters jHlavacek-Larrondo et al.| l |2012[ |. 
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Fig. 4.— Strong cool core (a = 2.75) bias and radio-loud AGN 
bias. The black, dotted lines represent the massiv e cool core 
“Phoenix” cluster, with radio spectral sl ope a.^ = 1.3 dMcDonald l 
let al.|2014|| and X-ray properties given inlMcDonald et al.| (|2Uil:i||, 
tor a range of redshifts. The pink shaded region represents a series 
of radio-loud AGN with spectral slope ag > 0, which demonstrates 
the range in possible AGN bias. The dark blue lines are the median 
bias in each bin and the light blue regions represent 2cr confidence 
intervals. In general, the cool core bias and the radio-loud AGN 
bias work to roughly cancel each other. 


[Strong CC 
Radio-loud 


- - - - Phoenix 


to our mock clusters (see §2.4): 


M; 


500 


\mJyy 


( 8 ) 


In Figure we show a contour plot of AC/C on 
Q^ooiM^oo, zj^a, Li ,4 parameter space. We display flux 
and luminosity data from a sam ple of X-ray s e lected 
massive elliptical galaxies from j Punn et al.| (| 2010 | 
and extreme AGN hosted in massive clusters lulavacekd 


(|2012|. As can be seen by Figure [Sj the 
e in SPT significance will be strong de- 


Larrondo et al. 

fractional change in SPT significance will be strongly i 
pendent on the assumed v alue of a.,, which h as a fairly 
large scatter. For example, Coble et al. (20071 measures 
a median value of a, = 0.72 between 1.4 and 28.5 GHz. 


Lin et al. (2009) explicitly measure the radio emission 
at trequencies ranging from 1.4 to 43 GHz in an X-ray 
selected sample of clusters and find mostly steep spectra 
(as > 0.5), but also find a substantial number of flat 
or even inverted spectra. Neither of these studies - or 
any published studies of radio galaxies in clusters, for 
that matter - extend to ~100GHz, so there is still sub¬ 
stantial uncertainty in the value of as that we should be 
using. 

In Figure]^ we show the strong cool core (a = 2.75) 
bias and the radio-loud AGN bias as a function of red- 
shift. Interestingly, while SZ surveys are biased towards 
cool cores, they are biased against radio-loud AGN, and 
this bias is of the same order of magnitude as the strong 
cool core bias, defined as the average A(//C over cool 
cores with a = 2.75. Since powerful radio galaxies tend 
to live in the center of strongly-cooling galaxy clusters, 
the net bias is partially canceled (see figure]^. In princi¬ 


ple, low-redshift galaxy clusters that harbor a powerful, 
radio galaxy (logLi ,4 > 25) but lack a cool core will be 
under-represented in SZ surveys. However in practice, 
such system s ap pear to be uncomm on in nature (see e.g. 


Sun| (|2009|) and |Dunn et al.| (|2010l)). The AGN bias is 
more stroiigly redshift dependent than the cool core bias 
due to the ^ I/di dimming where di is the lumi¬ 
nosity distance, whereas the evolution in the simulated 
cool core bias is a result of SPT beam effects and weak 
dependencies of our parameterizations on Pcr{z)- 


4. ROBUSTNESS OF SIMULATIONS 

In this section we estimate how relaxing several of our 
simplifying assumptions could change our bias estimates. 
Though a thorough treatment of each of these issues is 
beyond the scope of this paper (and should be addressed 
by, e.g. N-body simulations), an estimate on the relative 
importance of these effects can be obtained by recasting 
these effects into an M 500 dependence. 


4.1. Deviations from hydrostatic equilibrium and 
non-thermal pressure 

In the absence of complicated astrophysics, a 
galaxy cluster will equilibrate on dynamical timescales 
^dynamic ~ Rboo/v ^ ~ 10® yr. For hydrostatic 

equilibrium to generically hold, the cooling timescale 
should be much longer than the dynamical timescale 
tcooi ^ ^dynamic SO that hydrostatic equilibrium will 
be restored efficiently when gas cools. If the cooling 
timescale is shorter than the equilibrium timescale in the 
core of the cluster, we would expect pressure to be lower 
than what is required to maintain hydrostatic equilib¬ 
rium, since gas must be sinking towards the core. This 
would mean that the bias should be generically lower 
than what we have derived, since the observed change in 
pressure should be smaller than that required to main¬ 
tain hydrostatic equilibrium SPohs < SPeq. 

Throughout our discussion, we have assume that p res¬ 
sure is purely therm al, e.g. P = nkT. Analy tic work (Shi 
& Komatsu |20M I, numerical simulations ( Dolag et af. 


embolini et af.|2U13l [Nelson et af.|2Ui4[|Hatta^ia 

et al.|2014aD, and multi-wavelength observations ([Planck 
CollaDoration et af. 120131 [von der Linden et ah 


Donahue et ai. 2014 Sereno et al. 2014p, however, have 

shown that non-thermal pressure can contribute signifi¬ 
cantly (~ 20 %; increasing with cluster radius), whi ch can 
lead to an unde restimate of M 500 by about ^ 10% (Shi & 
Komatsu|2014 ). Sources of non-thermal pressure include 


turbulence from intracluster shocks, magnetic fields, and 
cosmic rays. If turbulence dominates the non-thermal 


pressure, non-thermal pressure will persist for on roughly 

the dynamical timescale (Shi & Komatsu 

2014). 

From equation ( 8 ) we see that bias sea 
to first order, we expect that adding non 

les as 1 /iW, so 
-thermal pres- 


sure will rescale the bias by (1 -I- P„on-therm/-Ptherm) \ 
if we assume that the processes responsible for forming 
cool cores and providing non-thermal pressure are un¬ 
correlated. As clusters evolve, this term could introduce 
an additional redshift dependence on the cool core bias, 
if the processes generating non-thermal pressure evolve 
with redshift. 


4.2. Varying fgas 
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Throughout this work, we have assumed that /gas = 
0.125 with no scatter. However, the numerical simula¬ 


tions of Battaglia et al. (2014b) which included radiative 
physics and AGIN feedback estimate a scatter ~ 10% in 
/gas at r < rsoQ. |Eckert et al.| (|2013|) suggest that cool 
core clusters and non-cool core clusters differ in /gas by a 
few percent, with cool core clusters more reliably tracing 
the cosmic baryon fraction. 

For a fixed number of baryons, a lower /gas will re¬ 
sult in an increase in the pressure normalization, in or¬ 
der to overcome the steeper potential due to additional 
dark matter. For a cluster in hydrostatic equilibrium, 
kT cx M/r giving P oc 

which means that the effects of changing l//gas will be 
roughly equivalent to changing M 500 . Equations (7) and 
(8) can therefore be rescaled by a factor (/gas/0.125). 
Adding intrinsic scatter in /gas would then simply add 
scatter in the SZ bias. If /gas varied with redshift, this 
could lead to a further evolutionary fa ctor in Equation 
(8). However, Battaglia et al. (2014b) found little evi¬ 
dence for such a redshitt de pendence, as has long been 
assumed ( White et al.||1993 ). 


4.3. Effects and Evolution of AGN 

We purposefully exclude any evolution in AGN prop¬ 
erties as a function of cluster mas s and redshift, despite 
evidence that such links exist (e.g., Ma et al. 2013). How¬ 
ever, by casting Equation ( 8 ) in terms ot the cluster mass, 
redshift, and radio luminosity, we allow the direct incor¬ 
poration of such trends into the bias estimate. By fully 
sampling a grid of parameters (see §2.4), we make certain 
that, regardless of how AGN evolve, we understand how 
they influence the SZ signal at all redshifts and cluster 
masses. 

Radio-mode AGN feedback can modify the gas distri¬ 
bution of galaxy clusters as gas is heated and expelled 
from the core, leading to deviations from hydrostatic 
equilibrium in the inner region. Such processes will typi¬ 
cally result in less pressure than predicted by hydrostatic 
equilibrium, since the gas is supported by feedback in ad¬ 
dition to thermal pressure. 

5. APPLICATION TO WELL-KNOWN SYSTEMS 

Using Equations 7 and 8 , we can now calculate how 
biased SZ surveys are (or are not) for some well-studied 
extreme systems. First, we consider two of the strongest 
kno wn cool cores: Abell 183 5 (7 = 0.85, M 500 ~ 10^® 
Mq; McNamara et al. 2006[ ) a nd the Phoenix cluste r 
(7 = 1.29, M 500 1.3 X 10 ^ Mq; McDonald et al.|2012). 

Assuming both of these clusters are at z = 0 . 6 , we hnd 
A/// =3.6% and 5.9%, for Abell 1835 and Phoenix, re¬ 
spectively. At z = 1.0, these biases are reduced to 2.8% 
and 4.2% (see also Figure]^. Given that these are two of 
the most extreme cool cores known, we expect the typi¬ 
cal bias towards selecting cool cores in SZ surveys to be 
<C 5% at z > 0.5. 

Figure quantifie s the radio bias fo r a variety of 
nearby clusters from Dunn et al. (2010) and Hlavacek- 
Larrondo et al. (2012), but here we specifically look at 
two well-known nearby radio-loud c entral galaxies: Hy- 
dra A (S '!.4 = 40.8 Jy, z = 0.055; iBirzaii et al.||2004|) 
and C ygnus A (S '!,4 = 1600 Jy; z = 0.056; Birzan et al. 
2004), the latter being one of the most powerful exam¬ 


ples of radio-mode AGN feedback known. At z = 0.6, 
these two clusters would be biased low in SZ significance 
by 1.8% and 110%, respectively. At z = 1.0, these biases 
are further reduced to 0.5% and 31%. Thus, while SZ 
surveys may miss the most extreme radio-loud clusters^ 
this bias is rapidly reduced with increasing redshift. In 
nearby X-ray-selected clusters, <Cl% of clusters have ra¬ 
dio luminosities as high as Gygnus A (Hogan et al. in 
prep), so we expect this bias to be small overall. 


6 . MASS BIAS IN SZ SURVEYS 

To translate SZ observations to cosmological con¬ 
straints, it is necessary to estimate mass distribution of 
an ens embl e of clusters. T o this end, [Vanderlinde et ah] 
(2010) and Benson et al. (2011) have adopted / as a 


proxy for M 500 by assuming a scaling relation / oc . 
As a consequence, the true M 500 of cool core clusters 
will be underestimated by an amount on the order of 
/ ^’(© 500 , 7 )(< 5 C/C)'^ 05 OO dy on average (for i? ~ 1 ) if the 
calibration was performed using only a sample of non¬ 
cool core clusters. In principle, this redshift-dependent 
bias should translate into a distorted N{M,z); however 
the log-normal intrinsic scatter in / given mass has been 
measured t o be 0.21 ± 0-10, c alibrated using X-ray ob¬ 
servations (Benson et al.||2013). Given that the Phoenix 
cluster has the most-cuspy X-ray surface brightness pro¬ 
file of the SPT clusters with Chandra X-ray follow-up 
(McDonald et al. 2013), which includes ~90 clusters, we 
expect the bias in / to represent an extreme example. 
Therefore, the scatter in / given mass due to cool cores 
should be a factor of several below the overall scatter 
in / and its current measurement uncertainty. The un¬ 
certainty in scatter has a negligible effect on the cosmo- 
logical constraints fr om current SZ cluster surveys (e.g., 
Benson et al. (2013|)), so the additional scatter from the 
effect ot cool cores will be even less significant. 

Although we have focused on the SPT SZ survey, our 
results should be applicable to other SZ surveys, e.g. 
Planck, AGT, etc. Broadly speaking, we expect less bias 
in a survey with lower angular resolution and greater bias 
in a higher-resolution survey, since convolving a wide sur¬ 
vey beam function with a cuspy SZ signal will smooth it 
to a less cuspy signal. 


7. GONGLUSION 

Using extensive Monte Garlo simulations of galaxy 
clusters and mock SPT observations, we have estimated 
the SZ bias due to cool cores in relaxed, massive sys¬ 
tems. By doing so, we have constrained the cosmolog¬ 
ical and astrophysical bias due to the presence of cool 
cores and radio-loud AGN in SZ surveys. We find that 
the bias from cool cores is no larger than ~ 10 % for 
z > 0.1 systems, and for typical high redshift objects 
(z > l,r 5 oo ^ 800kpc, Ay < 1) the bias is at the per¬ 
cent level. Further, the presence of radio-loud sources in 
cool cores should reduce the overall bias, though at low 
redshifts z < 0.3, the bias from radio-loud point sources 
should dominate any cool core bias. 

Our results support the long-asserted claim that an 
SZ-selected sample of galaxy clusters is a robust cosmo¬ 
logical probe: though we observe a small bias in the mass 
estimator, the magnitude of the bias is much smaller than 
the typical scatter in the mass relationship. 





















































We provide estimates of the SZ bias as a function of 
redshift, mass, and cuspiness, parameters that are model- 
independent, as well as the radio bias. One can estimate 
the bias of a given system easily by plugging in values to 
our function AC/C = f{z, rsoo, 7)+5(-^i.4, 2^, ^ 500 ) where 
f{z, r 500 , 7 ) is given b:^quation Q and 5 (^ 1 . 4 , 2 , M^qq) 
is given by equation By quantifying the cool core 
bias, astrophysical constraints on cool core properties 
from SZ surveys can now be more reliably interpreted. In 
particular, constraints on the cool core fraction from an 


SZ-selected dataset are only subject to a systematic bias 
of order one percent, a significant reduction over X-ray 
selection. Since there is a stringent upper limit on the 
redshift evolution of the cool core fraction bias, we can 
now confidently say that almost all observed evolution 
in the cool core fraction reflects genuine cool core evolu¬ 
tion. With the arrival of results from Planck, SPT, and 
the Atacama Cosmology Telescope and others, SZ sur¬ 
veys might be ideal for studying the mysterious balance 
between heating and cooling. 
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